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Abstract 

We present a systematic construction of integrable third order systems based on the 
couphng of an integrable second order equation and a Riccati equation. This approach 
is the extension of the Gambier method that led to the equation that bears his name. 
Our study is carried through for both continuous and discrete systems. In both cases 
the investigation is based on the study of the singularities of the system (the Painleve 
method for ODE's and the singularity confinement method for mappings). 



1. Introduction 



The investigation of the integrability of second order differential equations has been 
one of the most important enterprises in the history of integrable systems. Initiated by 
Painleve [1] and completed by Gambier [2] , it established the importance of singularity 
analysis as an integrability criterion. The results of the Painleve-Gambier investigations 
are of capital importance since they showed the existence of new transcendents, known 
since then under the name of Painleve. Overshadowed by this momentous discovery, the 
work of Gambier on linearizable systems did not receive the attention it deserved. The 
recent discovery of integrable discrete systems has led naturally to a critical examination 
of the work of the 19^*^ century masters. In particular, we have shown that it is possible 
to find discrete forms not only for the Painleve equations, but, in fact, for every single 
equation in the Painleve-Gambier list. The equation #XXVII of the list of 50 canonical 
equations [3] , which we decided to call the Gambier equation, was of course among them. 
Its discretization necessitated a thorough understanding of the Gambier approach. 

The key idea of Gambier (but we are conscious that the historical truth may be 
different) was to construct an integrable second order equation by suitably coupling 
two integrable first order ones. The latter were well-known: at first order the only 
integrable ordinary differential equations are either linear or of Riccati type. The Gam- 
bier equation is precisely the coupling of two Riccati in cascade (and it contains as a 
subcase the coupling involving one or even two linear equations). From the point of 
view of singularity analysis this coupling of two integrable equations is not harmless. 
Each of the equations has singlevalued movable singularities. However, the singularities 
induced on the second equation by the singularities of the solution of the first one (and 
which would thus look superficially as fixed) may lead to multivaluedness. This feature 
makes the application of singularity analysis mandatory. Its implementation leads to 
the (algebraically) integrable forms of the Gambier equation. 

In perfect analogy to the continuous case, we have introduced in [4] the Gambier 
mapping. The latter is a system of two coupled homographic mappings (which play the 
role of the discrete Riccati) in cascade. The integrable forms were obtained through the 
application of the discrete integrability criterion that we have proposed under the name 
of singularity confinement. 

In the present work we shall address the question of the construction of integrable 
third order systems in the spirit of Gambier. Namely we shall start with a second 
order integrable equation and couple it with a Riccati (or a linear) first order (also 
integrable) equation. This enterprise may easily assume staggering proportions. While 
at second order one had only two first order building blocks at one's disposal, at third 
order there are minimally 24 equations (the Gambier list) to be coupled to the two first 
order integrable ones. The situation is even more overwhelming in the discrete case 
since it is well-known that each continuous equation of the Gambier list may possess 
several discrete avatars. In order to limit the scope of our investigation we shall consider 
coupled systems where the dependent variable enters only in a polynomial way. This 



leads naturally to the coupling of a Painleve (P) I or II to a Riccati. 

Historically the coupling of a P equation with a Riccati was first considered by 
Chazy [5]. He examined an additive coupling of Pj with a Riccati. Starting with Pi in 
the form: 

w" = Gtu^ + z, 

he introduced a Riccati: 

y' = ay^ + I3y + \w + ^. (1.1) 

This coupling is additive as opposed to the one introduced by Gambler which is multi- 
plicative and assumes the form: 

y' = ay'^ + {P + Xw)y + 7. 

(In the case of the Gambler coupling w is the solution of a Riccati equation). Since the 
singularities of Pj are double poles [6/{z — ^o)^)? the only coupling that is compatible 
with integrability is the additive one. Assuming that a 7^ 0, we can put /3 = by a 
simple translation of y and Chazy found that the only cases where the leading singularity 
does not induce multivaluedness were when equation (1.1) assumed the form: 

1 - fc^ 

y' = ^^y^ + w + ^, (1.2) 

where k is an integer not multiple of 6. Thus 5 cases had to be examined as A; = 6m + n 
with n = 1, . . . , 5. Chazy found the following necessary integrability constraints: 



n = 2 7 = 

n ^ 3 7' = 

n = 4 7" = /U7^ + uz 

n — 5 7'" = A*77' + J^, 



where fi and u are specific numerical constants. It turns out that for = n they are 
also sufficient. For k = 6m + 1 the first condition appears at k = 7. In this case the 
constraint reads: 

2304 

7(5) = 4877'" + 1207'7" y-f^ - 24^7' - 487. 

5 

This equation has the P property and is thus expected to be integrable. Still it is 
interesting to point out that this equation is more difficult to solve than the one we 
started with which is of third order. 

Chazy offiers only a rapid comment concerning the case A; > 8. In fact, the con- 
straints obtained are necessary, but not sufficient for higher /c's. We have examined 
the first few cases beyond k = 7 using the same method as Chazy, namely singularity 
analysis (but unlikely Chazy our approach has profited from the existence of computer 



algebra tools). It turned out that none of the cases we examined satisfied the Painleve 
criterion. So, although this is not a proof in a strict sense, we can suppose that no 
integrable cases exist beyond the 5 identified by Chazy. 

2. Coupling of integrable second order ODE's with a Riccati 

As we have explained in the introduction we shall not attempt an exhaustive treat- 
ment of all 24 [2] (or 50 [3], or more [6]) second order equations of the Painleve/Gambier 
list with a Riccati. Instead we shall limit ourselves to the simplest case, namely equa- 
tions where the dependent variable enters in a polynomial way (instead of rational) . This 
limits the research to just three generic equations: Pj (already examined by Chazy), 
Pii and the linearisable, G5 (number 5 of the Gambler list), equation. Both Pi and G5 
have dominant singularities that are single poles i.e. w ~ l/{z — zq). Thus the ade- 
quate coupling is through a multiplicative Riccati. An additive coupling would lead to 
logarithmic singularities in the Riccati and thus to multivaluedness incompatible with 
integrability. 

2a Coupling Pu with a Riccati 

We start with the canonical form of Pu, namely: 

w" = 2w^ + ZW + fj, 

and consider the following multiplicative coupling: 

y' = ay^ + {nw + + 7. (2.2) 

A gauge transformation on y can be used in order to put /3 to zero. Next, we proceed 
to determine a, 7 through the application of singularity analysis so as to ensure the 
Painleve property for the system. Equation (2.1) has of course the P property, and the 
expansion of its solution around a singularity is: 

(J , ,4 

w= \- . . . + aAz — Zq) + 

Z- Zq 

where cr^ = 1 and 04 is a free parameter (the second one beside 2:0). The coupling of w 
with y must not lead to multivaluedness. Thus the coefficient n of the coupling must 
be an integer. This is only the first condition and, by far, not sufficient. In order to 
proceed further we expand y around the singularity zq and assume that y either has 
a pole at the same location or it is regular. Substituting the expansion form (2.3) we 
can compute the terms of the series of y and obtain the compatibility conditions for the 
absence of logarithmic terms in the expansion of y. We find thus the following condition 
for n = 1: 

7 = a = 0. (2.4) 



(2.1) 



(2.3) 



Let us point out that a = 7 = works for every value of n: w is simply related to the 
logarithmic derivative of y. We have furthermore, for n = 2: 



7' = ck' = 0. 



(2.5) 



For n — 3 we obtain 



a 



7 




+ az 



,3 




(2.6) 



Isolating 7 and integratinj 



once, we find: 



aa'" — 3a'a" — a'az — ko? = 



(2.7) 



and putting = ^ we find 



(2.8) 



Thus the logarithmic derivative of a satisfies precisely Pn (with a free constant k). For 
n = 4 we find a more complicated condition: 



Putting = (/}' we can integrate (2.9b) (multiplied by a) for the quantity 0^7 and 
thus obtain 7. Then (2.9a) gives a 6**^ order homogeneous equation for (f) and putting 
u = <p' / 4> leads to a 5*^^ order equation. This equation passes the P test and is thus 
presumably integrable but its integration is a more complicated task than the equation 
we started with, which is only of third order. 

For n = 5 we obtain again as a first condition a = 7 = which as we explained is 
sufficient. For n = 6 a first condition (as in the n = 2 case) is a' = 7' = 0. However a 
second condition appears. In fact, for n > 4 the free parameter of the expression (2.3) 
a4 starts appearing in the compatibility condition which must be identically satisfied. 
Thus for n = 6 we find the second the second condition either a = and ji = —7/6 or 
7 = and ji = 7/6. Thus this coupling works only for some particular case of Pn with 
a specific ji. For n > 7 we have not been able to find any integrable case, besides the 
trivial o; = 7 = one. In some cases it is even possible to prove the incompatibility of 
the constraints. We surmise that the multiplicative coupling of Pn with a Riccati does 
not possess any integrable case besides the ones listed above. 

2h Coupling the linearisable G5 with a Riccati 

The canonical form of the linearisable equation, G5 in the Gambler list is: 



3a77' + a'7^ + 3a(l + /i) + + 9^" = 0, 



(2.9a) 



3^aa' + io? + 37(1 - //) + 3oLZ + 9a'" = 0. 



(2.96) 



w" = —3w'w — + q{z){w' + w^). 



(2.10) 



The Cole-Hopf transformation w = —u'/u reduces (2.10) to a hnear equation u'" = 
q{z)u" . The function q{z) is completely free. Given this fact one can make two different 
couplings. The first is the 'standard' one where the solution of (2.10) for given q{z) is 
injected into a Riccati: 

y' = ay^ + nwy + ^. (2-11) 

The condition for Painleve property for n < turns out to be a = 0, while for n > 
it is 7 = 0. In both cases (2.10) becomes a linear equation (either for y or for 1/y) and 
the remaining free function (7 or a) does not produce multivaluedness. 

The second case of coupling is when q{z) is itself proportional to the solution of a 
Riccati. Thus the coupled system now becomes 

w' = aw'^ + I3w + 'y, (2.12a) 
y" = -3yy' - y^ + nw{y' + y^). (2.126) 

Only the case a 7^ needs to be considered since when a = equation (2.12a) is linear 
and thus does not have any singularities. Since a ^ we can take a = 1. As previously 
the coupling enters through nw with integer n since the singularity of (2.12a) is a simple 
pole. For n < the system has always the Painleve property and thus /? and 7 are free. 
On the contrary, forn > we have stringent integrability conditions. For n = 1, 2 there 
is no solution for 7 leading to the Painleve property for the system. For n = 3 we 
find as only solution /5 = 7 = 0. For n = 4 we obtain the condition: 

1=-^P'+P' (2.13a) 

and 

p" = 12pp' - 16/3^ (2.136). 

Putting j3 = —j^, (2.13b) reduces to 4>"' — and we thus have elementary expressions 
for (3 and 7. 

For n > 5 we can obtain the two compatibility conditions in the form of a higher 
order nonlinear system for /3, 7. It turns out that for the first few cases studied this 
system has the weak P property [7] . We have not tried to integrate these systems since 
their integration is more diflBcult than the problem we started with. 

3. Coupling of a second order mapping with a discrete Riccati equation 

Constructing integrable discrete systems in the same spirit as Gambler is quite 
straightforward once the basic ingredients are available. What is needed is a detailed 
knowledge of the forms of the equations to be coupled and a reliable integrability detec- 
tor. The second order mappings which play the role of the P equations in the discrete 
domain have been the object of numerous detailed studies and we are now in possesion 
of discrete forms of all the equations of the Painleve/ Gambler classification. The dis- 
crete integrability detector is based on the singularity confinement that we discovered 
in [8] and which has turned out to be of uptmost reliability. 



The coupling we are going to consider is a homographic mapping (discrete Riccati) 
for the variable y: 

_ ^ {ax + /3)y + {r]x + 6) 

(where y stands for i/n+i, y for yn and a, P, . . .,9 depend in general on n) the coefficients 
of which depend linearly on x, the solution of the discrete P I or II. (Wc shall not present 
here the coupling of the discrete analog of the linearisable equation to a Riccati. As 
a matter of fact this equation is the simplest non-trivial member of the hierarchy of 
projective Riccati systems, the discretisation of which was presented in full generality 
in [9]). The mapping (3.1) can be simplified and brought under canonical form through 
the application of homographic transformations on y. The generic form of the result is: 

^_ {ax + l3)y + l ^^2^ 



y + {^x + S) 

Non-generic cases do exist as well, and foremost among those is the linear relation: 

y{^x + 5)- y{ax + /?) - 1 = 0. (3.3) 

In what follows we shall examine in detail the coupling of (3.2) and (3.3) with either 
d-Pi or d-Pn (under various forms). 

How does one apply the singularity confinement criterion to a mapping such as 
(3.2) when x is given by some discrete equation like d-Pi or d-Pn? The first step consists 
in determining the singularities of (3.2). As we have explained in [10] the singularity 
manifests itself by the fact that y is independent of y. (We say in this case that y 
"forgets the initial condition" or "loses one degree of freedom"). The condition for y to 
be independent of y is just 

{^x + 6){ax + p) = 1. (3.4) 

This quadratic equation has two roots which we will denote by Xi, X2: they can be 
easily related to a, /9, 7, 5. The confinement condition is for y to recover the lost degree 
of freedom. This can be done if y assumes an indeterminate form 0/0. This means that 
x at this stage must again satisfy (3.4) and moreover be such that the denominator (or, 
equivalently, the numerator) vanishes. 

Let us assume now that for some n we have x^, = Xi. The confinement requirement 
is that k steps later Xn+k = ^2- (We must point out here that we require that Xn+k 
be equal to X2 and not to Xi again. If the latter were done this would mean that 
Xi is a singularity occuring periodically. Such singularities are not really movable, i.e. 
their position cannot be freely adjusted by choosing the appropriated initial conditions. 
Our conjecture is that they do not play any role in integrability, just like the fixed 
singularities in the continuous case). Starting from Xn = Xi and some initial datum 
Xn-i, we can iterate the mapping for x and obtain Xn+k as a complicated function of 
Xn-i and Xi. Since Xn+k depends on the free parameter Xn-i there is no hope for Xn+k 
to be equal to X2 if Xi is a generic point for the mapping of x. The only possibility 



is that both Xi and X2 be special values. What are the special values of this equation 
depends on its details, but clearly in the case of the discrete P's we shall examine here, 
these values can only be the ones related to the singularities. To be more specific, let 
us examine d-Pn: 

_ zx + a ,^ ^, 

X + X=- 2- ("J-^) 



The only special values of x are the ones related to the singularity Xn = ±1, Xn+i = 00, 
— T1 while . . ., Xn-2, Xn-i and Xn+3, Xn+4:, ■ ■ ■ are finite. This means that the 
two roots of (3.4) must be two of {+l,oo, —1} and moreover that confinement must 
occur in two steps. The precise implementation of singularity confinement requires that 
the denominator of (3.2) at n + 2 vanishes (and because of (3.4) this ensures that the 
numerator vanishes as well). Moreover, we must make sure that the lost degree of 
freedom (i.e. the dependence on y) is indeed recovered through the indeterminate form. 

3a Coupling various d-Pi 's with a discrete Riccati 

In this subsection we are going to analyse the coupling of four different forms of 
d-Pi to the homographic mapping (3.2) and to a linear equation (3.3). The d-Pi's we 
are going to consider are the following (presented below together with their singularity 
pattern) : 

x + x=^ + ^ {0,00,0}, (3.6) 

2: 

x + x + x = — + l {0, cxD, cxD, 0}, (3.7) 



x + x= — + 1 {0,00,1,00,0}, (3.8) 



X 



xx=-^-\-\ {g, 0, 00, 0, g), (3.9) 



x^ X 



with Zn = an-\-h and g^ = go A" (a, 6, go, A are constants). (More forms of the discrete 
Pi [11] are known but we shall restrict our analysis of the possible couplings to just 
these simplest forms). 

All the singularity patterns above have as common characteristic that one can enter 
the singularity through and exit it again through 0. This means that the condition 
(3.4) can have as a double pole. This results to the following conditions: 

ao + p7 = 0, ^ ' 

and since neither a nor 7 can vanish (lest the x^ term disappear) we have 5 =1/(3 and 
7 = One can, of course, consider the case where one (or two) of the roots of 

(3.4) are equal to 00: after all 00 is part of the special values of the singularity pattern. 
It has turned out that except for the case (3.8) the consideration of these cases does 
not lead to any interesting result. (Let us point out that the value 1 appearing in the 
singularity pattern of (3.8) should not be considered as a special value: it may well 



occur outside any singularity pattern). Thus the first discrete Riccati we are going to 
consider is of the form: 

_ y{ax + (3) + 1 

In all the cases considered, the first confinement condition, namely that y (at a suitable 
n) assumes the form 0/0 does not suffice in order to reintroduce the dependence in the 
initial conditions. It is thus necessary to proceed to the next order and introduce one 
further constraint (which turns out to be sufficient). Let us work out in detail the case 
of the d-Pi (3.6). Starting with a; = we obtain y = (3 i.e. independent of the value of 
y. For X = oo we obtain y = —(3(P and finally at the next step, x = 0, we ask that the 
numerator and denominator of y vanish. This leads to the first condition 

^13'^'^ =\. (3.12) 

Implementing this constraint leads to a second confinement condition that reads: a/ (3 = 
a/ p. This means that a = c(3 where c is a constant with an even-odd dependence. 
The solution of the constraint (3.12) is straightforward. Taking the logarithm of both 
members and calling h = log /3 we find the linear equation 

b + 2b + b = (3.13) 

with solution b = {p+ qn){—l)'^. Simple solutions to (3.12) can be obtained from this 
last solution. On the other hand just by inspection we can obtain solutions to (3.12) 
where /? is constant: /? = ±l,±z. 

The case of the "standard" d-Pj (3.7) can be treated along similar lines. The first 
confinement condition reads: 

PP^P^J=-1 (3.14) 



while the second becomes too complicated to be exactly solved. We prefer to proceed 
using one particular solution of (3.14) corresponding to constant /3's, for example j3 = i. 
This leads to a second confinement condition a/z — a/'z. Thus a = cz where c is a 
constant with ternary freedom (c = c). If we implement (3 = e^*'^/^ and define % = —az, 
we get, as a second confinement condition, the equation: 

X + X + X = 3/3'- + c, (3.15) 
X 

where c is a constant of integration. Thus, after considering the coupling with a d-Pj 
we get a d-Pj of the same type as one of the confinement conditions. This is in perfect 
parallel to the continuous case of Chazy (coupling (1.2) with n = 4) where we get 
another Pi as the integrability condition for a coupling between a Riccati and a Pj. 

The case (3.8) leads to still more complicated equations. One way to simplify 
them is to choose /3 satisfying: 

PI3'^P=1 (3.16) 



which is sufficient (but not necessary) to satisfy the first confinement condition. We 
can then implement the solutions j3 = i and /? = 1. If = z, the second confinement 
condition is a = c/2; (where c is a constant with quaternary freedom c = c). If /? = 1, 
we define x = —az and we get, as second confinement condition, the following equation: 

Az 

X + X = +c, (3.17) 

X 

where c is a constant with binary freedom. So, again here, we get a d-Pi of the same 
type as the one we started with as a confinement condition. 

For the case (3.8), it is also possible to consider a coupling where the condition 
(3.4) has and 00 as roots. This means that a = (we could also choose 7 = but 
these two cases are equivalent under the homographic transformation w ^ 1/w) and 
5 = 1/13. The first confinement condition then is 7 = —5. We define x — 1^1^ ^ind the 
second integr ability condition reads: 

X + X = -^ + l, (3.18) 
- X 

where c is a constant of integration. Thus again we get a d-Pi of the same type as the 
one we started with. Finally we can also consider the case where the condition (3.4) has 
00 as a double root. We then must have a = /? = 0. The first confinement condition is 
5 = —7 and we get the following relation for 7: 

(where A; is a constant of integration) which can be solved in an elementary way for 7. 

In the case of q-Pi (3.9) the full singularity pattern is one where we enter the 
singularity at q and exit it at q after four steps. However the complete study of this 
singularity pattern turns out to be untractable. Thus we shall limit ourselves here to 
the case where we enter the singularity through and exit it through after two steps. 
In this case the first confinement condition is just (3.12). Once this is implemented the 
second condition reads a/3 = Xa/3 which means a = cf3^^ where c is a constant with 
binary freedom and fx^ — X. 

Let us now turn to the case of the coupling of d-Pj with a linear equation (3.3). 
For the special values of d-Pj and 00, only three couplings have to be considered: 

y = ay+—, (3.20a) 
7a; 

y = axy+^, (3.206) 

y = . (3.20c) 

It turns out that in every case examined the second (3.20b) and third (3.20c) are always 
incompatible with confinement. The only remaining candidate is thus the coupling of 
the form (3.20a). By the appropriate gauge of y we can bring it to the form: 

y-y^—. (3.21) 
7a; 



Let us work out in detail the case of (3.6). A detailed analysis of the singularity pattern 
shows that if x vanishes like e then x diverges like ans x vanishes like — e. We 
compute the corresponding j/'s and find, at leading order, y ~ I/tc, y ~ 1/7^ and 
the condition for y to be finite is I/7 — 1/7 = i.e. 7 must be a constant with binary 
freedom (i.e. even-odd dependence). The analysis of the remaining cases proceeds along 
similar lines. For (3.7) we have the pattern {e,z/e, —z/e, —ez/z} and the condition for 
y to be finite is 7^ = 7^ i.e. 'j = k/z where /c is a constant with ternary freedom. 
The case (3.8) is related to the pattern {e,^/e, 1, —z/e,ze/z} leading to the confinement 
condition 72 = 72 i.e. 'j = k/z where /c is a constant with quaternary freedom. Finally 
the case (3.9) is related to the pattern {q + e, ae, — A/(a^e^), —ea/X,q} (where a is a free 
constant). Again we concentrate on the singularity induced by x = and which confines 
when x = again. This results to the condition 7 = A7 which means 7 = k/j,^ where k 
is a constant with binary freedom and iJ? — \. 

3b Coupling discrete Pu 's with a discrete Riccati 

In this section we shall examine the coupling of two different discrete forms of Pu 
with a Riccati: a difference one (which is the "standard" d-Pn) 

x + x=- ^, 3.22 

1 — X'' 

where z is linear in the discrete variable n and /x is a constant, and one of q'-type: 

x.= 4^, (3.23) 

x{x — 1) 

where q = qoX^ and p = poX^- 

Let us start with d-Pn (3.22). The singularity pattern of this equation is {±1, 00, =Fl}- 
This means that the singularity condition (3.4) must have ±1 as roots (the case when 
one root is 00 does not lead to anything interesting). As a result we have that d and 
7 are given by S — —j3/{a'^ — /9^), 7 = a/{a^ — (3"^). The pattern {1, 00, —1} leads to 
a confinement condition: 7 = {a + (3)a{a — (3) while the second pattern { — 1, 00, 1} 
leads to 7 = (a — l3)a{a + 13). Equating the two expressions for 7 we find gJ3 = 0.(3 i.e. 
j3 = ka where A; is a constant with binary freedom which we will ignore from now on. 
Expressing 7 in two possible ways we get finally for a the equation: 

This equation can be solved by linearisation just by taking the logarithm of both sides. 

The g-Pii has also two singularity patterns {g, 0, 00, 1} and {1, 00, 0, q}. Requiring 
and 1 to be roots of (3.4) gives the following expressions for 7, 5: 5 = 1/(3 and 
7 = • Next we obtain the confinement conditions for the two patterns of 

singularities: 

aa(3]3 + a(3_(3]3 + o^'^']3 + (3p'^']3-l = Q, (3.25a) 
aal3l3 + al30^ + ajBlB^ + 130^'^ -1 = 0. (3.256) 



Subtracting these two equations we obtain aj3 = aj3 i.e. j3 = ka where A; is a constant 
with binary freedom which we again ignore. Substituting back to (3.25) we obtain the 
final condition: 

which can be integrated through Hnearisation as explained above. The case where 
(3.3) has oo and q as roots is equivalent to the one treated above by a homographic 
transformation. 

We now consider the case where (3.4) has and q as roots which impose the 
relations S = 1/(3 and 7 = ^(g"*^^) • As first condition we then get that /3 is a constant 
with binary freedom. We ignore this freedom and consider /? as constant and we get 
the following relation for a: a = — -^^^g^- Finally, the case where (3.4) has q and 1 as 
roots has been studied but the resulting equations are far too complicated to be of any 
use. There is no other possible coupling of the form (3.2) with the q-Pn (3.23). 

Let us now turn to the case of a linear coupling given by equation (3.3). In 
the case of d-Pn (3.22) we require that the only singularities of the coupling terms 
{ax + (3)/{'^x + 5) be the two singularities ±1. This leads to a coupling of the form: 

y= xTi ' ^^^^^ 

where one of the parameters (e.g. 7) has been put to 1 through the appropriate gauge 
of y. Computing the successive y's we find that the condition for having a finite y, 
depending on the initial condition y, is just aa = 1. This means that all even cu's are 
constant while all odd ones are equal to the inverse of this constant. 

For q-Pii (3.23), in the case where (3.4) has and 1 as roots, we have two possible 
couplings: 

_ axy +1 (ooQ \ 

y = — j— (3.28a) 

and 

a{x-l)y + l ^ 

X 

It turns out that in both cases the confinement condition is the same as in the case of 
d-Pii namely aa = 1. When the roots are and g, the possible couplings are: 

axy + 1 /„ X 

y = — (3.29a) 

X — q 

and 

_^ aix-q)y + l ^ 

X 

The condition for intcgrability in the two cases is a = 1/A. Two other couplings are 
possible when the roots of (3.4) are q and 00: 



(3.30a) 



and 

y = a{x-q)y+l. (3.306) 

The integrability condition for (3.30a) is a — q and for (3.30b), it is a — 1/q. Finally if 
(3.4) has 1 and q as roots, the possible couplings are: 

- Q!(a; - l)y+ 1 

y = (3.31a) 

X — q 

and 

_^a{x-q)y+l^ (3.316) 
X — 1 

In the case of (3.31a), the integrability condition is 

aaa{—q + 1) + a{—qq + q + q — l) + q^ — q — (3.32) 
and in the case of (3.31b), the condition reads 

aaa{q^ — q) + aa(—qq + q + q — 1) — q+l = 0. (3.33) 

Equations (3.32) and (3.33) are integrable and they belong to the family of linearisable 
equations [10]. 

One last remark is necessary at this point since we have seen that almost all the 
equations we obtained contain terms with binary, ternary or quaternary freedom. The 
presence of these terms indicates that our systems must be augmented by adding more 
components. This will not alter the order of the resulting equation: it just increases the 
number of its parameters. The continuous limit is, of course, affected by this choice. 



4. Conclusion 

In this work we have presented a systematic approach for the construction of 
integrable third order systems through the coupling of a second order equation to a 
Riccati or a linear first order equation. Thus we have extended the Gambler approach 
(first used in his derivation of the second order ODE that bears his name) to higher 
order systems. We have applied this coupling method to both continuous and discrete 
systems (given that we have already presented [4] the discrete equivalent of the Gambler 
equation). 

One point remains to be discussed. It is often argued that, since the Riccati is a 
linearisable equation, the coupling of the Riccati to another of the same kind or to an 
integrable second order is always integrable. The (naive) argument is the following: first 
solve the first equation, substitute the solution into the second and solve it by linearising 
it. The argument about singularities is usually swept aside by the statement that one 
is interested only in solutions on the real-time axis. However the situation is not that 
simple. What integrability consists in is a global description of the solutions of the 
equations. The argument about solutions on the real-time axis is not acceptable since 



it offers just a local description of the solution of the equation. A global representation 
of the solution of a linear equation (and, thus, also of a Riccati) involves path integrals 
winding over the complex-time plane. Thus the study of movable singularities is crucial 
and the Painlcve property a necessary condition for integrability of the systems. 

How do these arguments carry over to the discrete setting? One must go back 
to the way difference equations are formally solved. Given a linear difference (or q—) 
equation, we can express the solution as an infinite product of matrices, the elements of 
which depend on the coefficients of the equation. A singularity appears whenever one 
of the matrices is singular. In this case the solution of the linear difference equation 
cannot be defined for every n. However it is in general possible to choose the coeffi- 
cients of the equation so as to avoid these singularities. In the case of a coupling the 
coefficients depend on the solutions of some other equation. Thus there is no way to 
control the singularities (which depend on the initial conditions of the first equation). 
As a consequence the solution of the second equation cannot be defined everyTvhere 
unless the confinement property is satisfied. Thus, again, despite the linearisability of 
the discrete Riccati, whenever we talk about a global description of the solution of the 
coupled system, the application of the adequate integrability criterion is mandatory. 
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